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ON TRANSONIC FLOW PAST A WAVE-SHAPED WALL! 


By CARL KAPLAN 


SUMMARY 


The present report is an extension of a previous investigation 
(described in NACA Rep. 1069) concerned with the solution 
of the nonlinear differential equation for transonic flow past 
a wavy wall. In the present work several new notions are 
introduced. which permit the solution of the recursion formulas 
arising from the method of integration in series. In addition, 
a novel numerical test of convergence, applied to the power 
series (in the transonic similarity parameter) representing the 
local Mach number distribution at the boundary, indicates 
that smooth symmetrical potential flow past the wavy wall is 
no longer possible once the critical value of the stream Mach 
number has been exceeded. 


INTRODUCTION 


In NACA Report 1069 (ref. 1) the writer considered the 
problem of two-dimensional transonic flow past an infinitely 
long sinusoidal solid boundary. The problem was treated 
in the physical plane and the purpose was to investigate 
whether or not the flow past the wavy wall remains a smooth 
symmetrical type of potential flow when the undisturbed- 
stream Mach number exceeds its -critical value. By a 
smooth type of potential flow is meant one for which the 
velocity potential, say, and its first derivatives are single- 
valued and continuous; that is, there are no discontinuities 
of the nature of shock waves. 

Initially, the Prandtl-Busemann small-perturbation method 
‘ was applied and the velocity potential developed inclusive 
of the third power in the “thickness” coefficient a 
where 
a amplitude of wavy wall 
wave length of wavy wall 
The velocity potential was then referred to the critical speed 


== 3/2 
of sound ce, the coefficient e was replaced by ka — My 


y 

h pa Ite | is th ic similari d 
where k = OMY” is the transonic simuarity parameter, an 
terms involving powers of 1—M „° higher than the first were 
neglected (main assumption for transonic flow). This simpli- 
fied or transonic form of the Prandtl-Busemann solution was 
shown to be identical with the one obtained by solving the sim- 
plified nonlinear differential equation for transonic flow, with 
the boundary condition taken not at the wave-shaped wall but 
at the flat plate corresponding to vanishing amplitude. The 
calculation was carried through the sixth power in the tran- 


sonic similarity parameter k and corresponds to the in- 
superable task of obtaining the Prandtl-Busemann solution 
to the sixth power in the thickness coefficient e. Thus each 
iteration step of the Prandtl-Busemann method contributes 
to the transonic form of the solution, which may therefore be 
considered a result of thin-profile theory with disturbances 
not necessarily small compared with 1—M,?. The main 
conclusion reached in reference 1 was that the transonic 
similarity parameter k must be less than <a value still 
somewhat greater than the critical value 0.83770 calculated 
there. 

The purpose of the present work is to express the solution 
of the problem of transonic flow past the wavy wall in a 
form more suitable for general considerations and to prove 
that the assumed smooth symmetrical type of potential flow 
cannot exist at stream Mach numbers beyond the critical 
value. 


SYMBOLS 
L,Y nondimensional rectangular Cartesian coordinates 
a amplitude of wavy wall 
A wave length of wavy wall 
Af2v reference element of length 
€ "thickness" coefficient of wavy wall, ax 72 
Y ratio of.specific heats at constant pressure and con- 
stant volume 
U stream velocity 
M local Mach number 
M, undisturbed-stream Mach number 
«tenons (ytDe 
k transonic similarity parameter, GMF” 
C;  -critical speed of sound 
$ velocity potential of flow 
f(x,y) transonic disturbance potential 
Fly) functions of y only, related fo f(z,y) 
As? numerical coefficients l 
A”? generating functions of k, SA? E a 
i T0 
Axes functions of dummy variable r, > n>? 7", lower 


n=1,m 
label starts from 1 when m is negative and from 
m when m is positive 
Primes denote differentiation with respect to independent 
variable. 


1 Supersedes NACA TN 2748, “On Transonic Flow Past a Wavo-Shaped Wall” by Carl Kaplan, 1952, 
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ANALYSIS 
GENERAL FORMULAS 
Tf the undisturbed stream is in the direction of the positive 


z-axis, then the velocity potential $ referred to the critical 
velocity Ce can be written as (see ref. 1) 


1 
y+1 


where the second term on the right-hand side is a disturbance 
velocity potential and implies that terms involving powers of 
1—M,,” higher than the first have been neglected. The differ- 
ential equation for transonic flow satisfied by the function 
f(z,y) 1s obtained from the general differential equation for 
compressible flow and takes the following simplified non- 


linear form: 
of of o a) 
Oy! Oz Oz? 


The boundary conditions to be fulfilled by f(x,y) are as 
follows: 


geet el — M ,?)f(z,y) 


. (at y= œ) 


of. 
ay 9 (2) 


Of : 
ov sin zx 


Here x and y are nondimensional rectangular Cartesian 
coordinates simply related to the physical plane coordinates 
X and Y by means of the transformation 


1—X 
y-(—M, Y 


^ - 


Ree ere ra 


and the equation of the wavy wall is Y=e cos X or 
y=(1—M,,”)"*e cos x. Clearly, f(x,y) involves only the 
variables x and y and the transonic similarity parameter k. 


The most general form for the function f(z,y) to ensure 


symmetrical potential flow past the wavy wall is the following 
(see fig. 1): 


fG3)——24 21. sin nz (3) 
where the f, are functions of y only. When this form for 


f{(z,y) is substituted into equation (1) and the coefficients of 
the separate harmonic terms sin nz are equated to zero, 


Y 
X f PAANAN = c NN Ah : z AEEA E X 
Y 


Figure 1.— Wave-shaped wall. 


Y 
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the following system of second-order nonlinear ordinary 
differential equations for f, results: 


f, —mfe- nim (n—m) fuf m — 
g^2jm (nM) Fu fam (n=1,2,...0) (4) 


Before proceeding to the solution of these equations, several 
formulas of general interest and subsequent use are given. 


They have been derived in reference 1. "The local Mach 
number distribution is given by 
of ; 
— M?= — (1 — = 
1—M?-—-—(1— M. 22 (5) 


The equation from which the critical value of the transonic 
similarity parameter is calculated follows from equation (5) 


by taking z—0, y—0, and M=1; that is, 


OfN — 
(35 "A 


or with the aid of equation (3), 


25nf,(0)—1 (6) 
"The pressure coefficient 
C, M wo P— Pe 
| g Pa” 
is given by 
m 2 of 
C, «mu 07 MAE) (7) 


INTEGRATION IN SERIES FOR THE FUNCTIONS fn 


In reference 1, equations (4) were solved by an iteration 
procedure and approximate expressions for f, to fe were ob- 
tained. An examination of these expressions showed that 
the general form of f, is 


iia. 
fepe cem SS yt Sy Ar? kt (n—1,2,...9) (8) 
p=0 g=0 r=p 


where, if y=0, the upper limit of q is n—1 and, if p +0, the 
upper limit of q is 2p-++-n—2. The four-labeled coefficients 
A; ?are real numbers calculated from recursion formulas 
obtained from the system of differential equations (4) and 
the boundary condition at the surface of the wavy wall. 
The boundary condition at y= œ is automatically satisfied 
by the form of f,; whereas the boundary condition at the 
wall t&kes the form 


wr m (6) 


Inserting the expression for f, given by equation (8) into 
equations (9) yields immediately the following results: 


09—1 
2 (2p-+n) Az P 5 Az? (n=1,2,... K (10) 
p=0 2501 


ON TRANSONIC FLOW PAST A WAVE-SHAPED WALL 983 


where, if n=1, the lower limit of p on the right-hand side 
is unity and, ifn +1, the lower limit of piszero. Also, if n=1, 
the upper limit r of p goes from 1 to œ and, if n1, r goes 
from 0 to o. 

By elementary manipulation of series, the second of equa- 
tions (10) can be replaced by the following more useful forms: 


nAz $= Ar? (n=2, 3, . . . œ) 
i. Q (r=1,2,...0\¢ OD 
: PE. !— 2:0 p--248: fie. ot 


where in the first term on the right-hand side the lower 
label p starts from 1 when n=1 and from 0 when nz. 

In reference 1, recursion formulas were derived for the 
coefficients A"?. In the present report a much more 
significant approach is introduced. Note that equation (8) 
can be rewritten in the following form: 


n—l 
fa E lecto SA 
p=0 g=0 
where . 2 
A^? m 214,4 ? ANE 


“AN? y* (12) 


In a manner similar to that described in reference 1 for the 
coefficients A7 ?, recursion formulas can be obtained for the 
power series Á”. Indeed, the two types of recursion formu- 
las are intimately connected and, for a given value of y, 
the one can be obtained from the other by mere inspection. 
A single recursion formula can be written for the general 
quantity Á”, but the resulting expression is cumbersome 
and serves no practical purpose. It is much more desirable 
to obtain separate recursion formulas for each value of p. 
As examples of procedure, the recursion formulas for p—0, 
1, and 2 are considered in the sections which follow. 


"RECURSION FORMULA FOR Anr, 


With p=0, the recursion formula is (compare with eq. (57) 
of ref. 1) 


2n(q+ Az 2=(g +1) (g-+2)8" z ARs 
eig (m4-1) (n—m—1)4*-7*- sei. 


2 d anra 


(n—2,3,...9; q—0,1,.. =. (13) 
where 
0 (q—n— 2) 


1 (q £n — 2) 
This recursion formula can be solved, the solution starting 


with g=n—2 and descending towards Lied Thus, for 
g=n—2, equation (13) becomes 


ô” c mes 


8(n— 1).47 4$ (m+ 1)(n—m—1)A" o An n1 
(n=2,3,...0) (14) 


Now, multiply both sides of this equation by 7* where 7 is a 
dummy variable and sum from n=2 to n— o. Then 


o * ceo n—2 
82 Yn— 1) A; tr^ —2 jm 2 (m+) (n—m Am An nct 
(15) 


Let l 
A, 2 Ar 
then 
Ay =D nA? 2 77i 
and = 


(d n—1l 
(4 oP Dy 77 $3 m(n—m) Az. 1 Anon t 
By observation, it can be seen that the right-hand side of 
equation (15) is equal to (rA, ,?. Equation (15) can thus be 
replaced by the following first-order nonlinear differential 
equation: 


(rA, o) —8rA, o +8A,; ,—0 (16) 


The solution of this equation is 


A, ) —c,e* qr o 
where c, is the arbitrary constant of integration, From the 
definition of A; o it follows that, with r—0, 

Co= Ai,” 

and hence 
5 1l, gar? 

Ayo a^a 

ae a 


Note that by definition A*,° involves the coefficients Az? 
which are ultimately calculated by means of the auxiliary 
relations, equations (11), engendered by the boundary con- 
ditions at the surface of the wavy wall. In solutions of the 
recursion formulas, therefore, the coefficients A*,° appear as 
undetermined quantities. 


The expansion of zs > in powers of rA! is a nonirivial 


problem which fortunately can be solved with the aid of 
Lagrange’s investigations on the reversion of power series. 
In reference 2, Lagrange’s problem is illustrated by the 
following example: 

To expand e* in powers of y=ze"., 
as 


The result is expressed 


u TNT 
= 1 ray 450 20) En) jig ey m 2) JP. xr 
a(a —nby 


"| 


y - Bes (18) 


In particular, with a—1, b— —1, and £—e*, the solution of . 


- the equation log {=y or =e (generally referred to as 


Eisenstein's problem) is given by the series 


= -5 


Then, if r is replaced by 
of equation (17) is 


E ç and y by qr Ate, » the solution 


A E o n^ 3 
Ay apie) 


n-i p” -1 
n=l í 


984. 


or 


24 AN i um 


Hence, by equating coefficients of equal powers of r on both 
sides of this equation, 


Zuge. ort 


‘At §=— = (At, ye (m=2,3,... æ) (19) 


nid" 
which is the solution of the recursion formula, equation (14). 
The coefficient A*_{ can be considered as the. generating 
function for the set of coefficients 4,*,°. Thus, equation 
(19) can be written as 


eo n* —2 ro) a 
x O7£,2 l-. 10 

X AS (2 are) (20) 

Then, by equating coefficients of equal powers of k on both 

sides of this equation and putting Aj$=1 (see boundary 


condition, eqs. (10)), the following equations are obtained: 


n* 2 
Ant, Q n|4* -! 
ut NN. dE A19 
D1— (n—1)4^7! - 
jc 10 NE Su " 
- En 12 (n—1)4*^7! 4 rta 2)4*-1 (Ao i? 
Note that from the first of these formulas 42 i- and that 


from the boundary condition, equations (11), A? =r 


Consider now g=n—3. Then equation (13) becomes 
2n (n— 2) AS -5—(n— 1) (n—2) Az -1 


sum (n—m) Az _$ Asin! (22) 


Multiply both sides of this equation by r* and sum from 
n=3 to n—o. Then 


2 2; n(n—2) Az -gr*—5 (n—1) (n—2) Az 37^ 4- 


e n—l 
321 nr $im(n—m)Az-$4z-z.1 Q3) 
Let 
A, =D) Arr 
and 


A, = 2) A: zar 


REPORT 1149—-NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


then 


^ End 0 


Ayo Az =>) r2 m(n—m-+ 2) Àz.1 


It can be shown easily that the last term on the right-hand 
side of equation (23) is equal to 5 r (r14 o Aa dY. 


equation (23) can be rewritten as the following differential 


Hence, 


- equation: 
277 Ag y — 2r Ás y =r A, y — 2rÀ; 43-24; ter (Ay p Aso)! 
(24) 
Now, from equation (16), it follows that 
(rA. o =8 (ri o — A) o) 
_ 4A d 
(rA; o) 4—r4, ; (25) 
rA, o (r. A, o) —4r4À; 0” 
Then equation (24) can be written as 
; 8+r(rAi 6) NN 1 N 
(As a pA A 0429) 7 Rara, gh) A) 
6rAı o +441 c] (26) 
Let rA: o =v m where v is the new dependent var- 
iable. Then equation (26) becomes 
= =i r(rAi Jy [2r (rAy o)! — 6r Á, o t4A, o] 


or with the aid of equations (25), 


v’ =— 5 At o) 
Hence 
v ps T Ay o’ 4l c 
8 10 8 1 
or : 
(r.A; o» EU TT 9). 
8r dro = g m TA, o ry n Arii d o 


The arbitrary constant c; is determined by differentiating 
this last equation twice and evaluating for r=0. "Thus, 


A? 0 
$64 ary 
Finally, 


Aa o an [(Ay 0) 4-8 -rra e s TAY v. 0 (27) 


Now, from EURO (18) with a=2, b—-—1, and f=e?, it 
follows that 
n~2 
—22m—2] a IY 


n-2 
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Or with ? 
: I 
[= =" To mr rA’, 
eo xz —3 
(Ar 223 uei Ayr 
and 


e qm -3 _ 
(4: V= 2 > gji (AN) r7 7? 
Thus, equation (27) becomes 


2jn4r.i E. p^-7lol =5 P2 Ge ns -5 (A! oe pice 


ES Li n 0,n2—1 
8 aly (n DAR -37* 
Hence, replacing 43 1 by its value from equation (19) and 


equating the coefficients of equal powers of r on both sides 
of this equation yields 


—2 20 "n^ 
Ao 5 (n— a= (At) 


(n—3,4,...c) (28) 
Again, 


A; -2—8 zi CAT 


(n— Oe 


the solution of the recursion formula, equation (22). 
equation (28) can be written as 


3A rea LS e 214 Ai? * (3:4it p") ey 


Su mes (lae (29) 
2(n—3)4^ NES ?* 
Then, equating the coefficients of equal powers of Æ on both 
sides of this equation and putting 44 3=1 and Aj fa 
yields: 
, n* -3 n* -4 
Taice a a 
Anta ae 1 aaa 
8n"? m3 
A,tal=G oop 438+ Gogg At 
8n" -3 ` 
Aane naaa T (30) 
a sina 
(n—3)4^71| e? e 


AQ P4 | ss aes (Ag D? 


Note that from the first of these formulas A? 9 T and that 


from the boundary conditions, equations (11), 43 0—7y 


The procedure followed in order to obtain equation (28) 
for 42.3 is a general one and with very little difficulty other 
members of the family A can be obtained. For example, 


Az -$—48r ayaa ead 99)" 343,9 4- 


1 n -—4 ‘ 
5 CDT AU AHANY] (81) 


From this equation it follows that 


" "ONE net 1 mnt 
Ay -3 97 (n—3y4*3 t3 (n—4)47-! 
A," ogg o A3 9-8 BEEN A204. 
30177 (n—3)4*717 9?! (n— 444-1 1 
(32) 
a ae Eo 
4*17?!| 3(n—4) 2 (n—5)! 
From the first of these ra and he baid conditions, 
ations (11) í 
rd » AL o— ii - 1536 


At this point, it is noted that the coefficients of the form 

* 0 are calculated from the first formula of each set, equa- 
tions (21), (30), (32), and so forth. A number of this type of 
coefficient have been evaluated (see ref. 1) by means of the 
recursion formulas for the coefficients A? ? themselves. 
They are listed as follows: 


4i =a di =i 
Ai =e ig UEM 


À. careful examination of meee numerical values leads to the 
general rule, 
3n— 5 
At 8 =i (n=2,3,... œ) 
and from the boundary conditions, equations (11), 


{3n—5} 


ani (n—-2,8,... e) 


E 
where by definition 
(39n—5]—1X4X7X10X13X ... X(3n—5) 


In the expression for the local Mach number distribution 


936 


evaluated at the crest of the wavy wall (z= 
occurs the following power series: 


0, y—0), there 


F=>\nAg 3k 
n=l 
or 


q4^ 71 


Fei RE {3n—5} za 


This power series can be expressed in the closed form 
P=2 Sg) 
-2-:-ii) 


The graph of F against k is a semicubical parabola with the 
cusp point at k=% and F=2. With the necessary condition 


that one and only one value of k correspond to & given value 

of F, the transonic similarity parameter k cannot be greater 

than 4/3. Moreover, the lower limit of k is zero when the 

amplitude of the wavy wall is zero but M, is different from 

unity. ! a 
RECURSION FORMULA FOR A »4 


With p=1, the recursion formula for 4^, is (compare 
with eq. (69) of ref. 1) l 


4(n-- 1)4"7 — 2(n--2)(a HA Hla HLH 2A 
-3ni By OSD namang Ag cnn Asta 


2 Q,70 m»-g— 4, 


(n—1,2,...9; q—0,1,...-) (33) 
where 

(qq=n—1 or n) 

l (q(=n—1 or m) 


. 0 (q—7) 
a1 (qxn) 


The solution of this recursion formula proceeds as in the 
case p=0, the solution starting with g—7 and descending 
towards g=0. Thus, for y=n, equation (33) becomes - 


Atria -g^ AX ATH? — (n—1,2,...0) 


or 


AMAS——IG—1)4M Ar} (n—2,8,...9) 


Hence, inserting the expression for 42.9 given by equation 
(19) gives 


n-3 i : 
= (A)**  (n-2,3,... 


Ae $m 2)l4" 


c) (34) 


or 


Í q^ -3 


a M O Nati 
24 Ali doi e (2 ite) 


emn e ————— € ——— . 
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From this relation, by equating coefficients of equal powers 
of k on both sides of the equation, the following equations 
are obtained:  ' 


1 qm -3 
Aii =g Ge 
n” -3 1 nm -3 
2&—11.. 10 —— Ale 
Aiii 2 maipi 40-3 Gaye 4 


T -— 
: GE [45 $-F CA3 2] — (85) 


1 n^? 101. of 
8 (n —3)i47i [Aj 8+3(4A5 2] — 


Asai I= — 


TR (n=2,3,... c) 
Gondar now g=n—1; equation (33) becomes 
4(n-+ 1)Àz -i= 2n(n t 2)47,' — 


3 (m--1)(n—m —1)A**5 Araia 
5 n(n-F 1439 uu (n—1,2,... œ) 


Multiply both sides of this equation by r* and sum from 
n=1 to n=æ. Then 


- 43) (n+1)Arir= 22 yn(n-- 2)4*,! r° — 


32 (m1) (n—m—1) Amt * Anc tha 


= = AY y nint 1)4^, r^ (30) 


=>) Ar 1 
^ fa] 


= Dy Az_or A, =>) An gr 


Then it can easily be shown that the second term on the 
right-hand side of equation (36) is -$ T Ao ,À, ,) and 
that equation (36) can be replaced by 


4(rA, y= 2r* A, vt 6r A, /—5r( A, JA, ne râ 10 Ay d 


37 
Now, oe 


Ag 1=—= At? $5 (n— 1)42 2 17 
' 8 fn] 
or, with the aid of equations (25), 


T Á, =- Al, (rÁ, yy 
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Pd 


Hence, from the use of equations (25) and equation (27) for Aso’ it follows, after routine calculations, that the integral 
of equation (37) is 


rÅ = Ay i-i 3 (rA -I Z (r4, o (rá ECCE) E (7 Ay ,) — 2r* A, o” I 


where 
Ar 
01—64 7710 (4^, oja 
Now from equation (18), 
^m Smet R , 
(rA; o P= 25 v mim Ve ro (38) 
Then ; 
A Ar- pml (A1 P n7? peni 1.2 nr -5 
2j 43 regag A "1-23 gue" ideo bape” t 
nr 7? 
+052 gone”? Beene" | (09) 


Equating coefficients of equal powers of r on both sides of this equation gives 


Ant — READ Glo got e el egre n apte ouis Le T gs oats 


(n=2,3,... e) (40) 


From this equation the expressions for all the coefficients of the type .42-i.1, with n=2, 3,... e and r=0, 1, 
. œ can be obtained. Thus, for example, i 


1 ne -5 


Anna = — 32 a LN gte pitu i] 
Arni Cere At a "P (ome DA He3240 Boy gr 7-8 A + 


. (41) 
256 A3] c ep A 14-848 9c uy 
For gq=n—2, equation (33) becomes 
4(n-+1) Ai = 2 (n— 1) (n - 2) Az _i—n (n— 1) A" i JS [(m 4- 2) site 2) Amt? 0 4»---21 
(m-+1)(n—m—1) Ant * Anni n(n) AP AnH? (42) 


The solution of this recursion formula follows along the same lines as that for equation. (36) and leads to the following result: 


1 | 102 414 205 9 1 ns 


A = Fie n St Goat acu Hert) | ar etiem By gaat 


12128 (is gr 2 ap pats eH arte e] m 
" (n=38, 4, . œ) (43) 


From this pee there follow, in the usual manner, expressions for coefficients of the type Anzio -+1 Thus, for 
example, with r=0 


aml li. S E n? 3 AL ; n? a 8 205 EC M m7 -7 


T (44) 


988 
RECURSION FORMULA FOR 4*4 


With p=2, the recursion formula for A”, is 


8 (n+ 2)A* ?— 2(n4- 4) (q+ De 1? AR 
(q3-1) (q-T-2)57 *17 *? Ans 


"— "SS" m —2)(n—m4- 


q1::0 Timg—gi 


EL tats D» $4 m(n— 


m)" ani n(n--1) (32 AM An 1 

pn he? Alt Antio 50.042 411 Artt o) 

n(n--2) ^3 AA? An *2 9-38 A29 An 22)) 
(n=1,2,.. 


.0; g=0,1,...n+2) (45) 


where 6 is defined in the usual manner. The solution of 
‘this recursion formula proceeds as in the previous cases, 
starting with gq=n-+2 and descending towards q—0. Thus. 
for gq=n-+2, equation (45) becomes: 


Ani——i nA Anse (n=1, 2,. ui .«) 


or 


Anim (n—2)439 A22 — (n—83,4,...c) 


Then, from equation (19), 


Ana -za Ne aF — (A},0)* +2 
^ (n=3,4,...0) (46) 
or i 
co (n—2 212 
AnA k= on (as: ey 


By equating coefficients of equal powers of k on both sides 
of this equation, the following equations are obtained: 


n—2)n* 7? 
ajm en 
n22 (RM 4) 7 4, 
ae (47) 
2— 4)n* - 1 
Anzi 0 EE “| aite aA | 
Consider now g=n-+1; equation (45) becomes 
8(n4-2) Az ,1— 2(n-- 2)(n-- 4)47 ,$— 
n—2 s 
-g^23(n-- D —m— 1) Ami Anny 
z nn HINARANA AnH 9)— 
nint JAPA + A?,°A* 12 0) 
(n—1,2,... Œœ) (48) 
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With the introduction of 


Ao 
A a= 2 4; +1 T” Áo = 254"! r 


equation (48) is replaced by the following ordinary linear 
differential equation: 


S(r A a)’ — 2r*(r Á s P) -12r(rÀ 4 ,)2- 10704 3 ,— 


FPA, o s 5 Ayr? A, 5 APA, c — 


A99 rA: s" 4-49, A, o —A7 74,4 S" + 

A? As pu (49) 
By repeated use of equations (25), the solution of this differ- 
ential equation is found to be as follows: 


8r*A , ,— (A! vis z CA )— dd (rA, Pa (ra oP — 


aE (rA oY i E árg A, ot 


Bq 69 [gea vna of Irene 


* 


An | 4. EA Y Iran (50) 


Finally, with the aid of equation (38) and the definition of 
4a 2) 


Anz? ings ien" (n— R2 wm 
1 n^ -5 
^ BGR iei gis 69 [enr 
1 n7 1 OYn +2 
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A?,° 
(n=3, 4, 2... 0j €1— 64 ru) (51) 


From this equation, there follows*in the usual manner 
formulas for the coefficients of the type Az-2,3,. Thus, 
for r=0, ~ 


Az-i3 m3 qm Ts 
Agri a= =| ea -p pite a 3yi | 4599 (n—4)ld à * 
1 qm 
aaa Fe a a 


DISCUSSION OF CONVERGENCE OF SMOOTH TYPE OF 
POTENTIAL FLOW PAST WAVY WALL 

In the preceding sections, a number of examples of recur- 

sion formulas and their solutions have been given in consider- 

able detail. The purpose of this exposition is threefold: 


ON TRANSONIC FLOW PAST A WAVE-SHAPED WALL ; 


First, to show the inherent elegance of the method of integra- 
tion in series although the equations concerned are nonlinear 
in character; second, to present & type of analysis which may 
be useful in other problems involving nonlinear differential 
equations; and third, to indicate that an analytical proof of 
convergence may ultimately be obtained by careful examina- 
tion of the recursion formulas for the quantities A*,? and 
their solutions. One example is the obtaining of the general 
expression for the coefficients 78 and the subsequent con- 


clusion that k< = 


In actual practice it has been found niore convenient to 
evaluate the coefficients 47? from their separate recursion 
formulas rather than to derixe the general formulas. The 
appendix contains the exact numerical values of the coeffi- 
cients necessary for the development of the functions f, to 
the eighth power in the transonic similarity parameter k. 
These coefficients are utilized to demonstrate numerically 
the test for covergence of the smooth symmetrical. type of 
potential flow (eq. (3)) assumed in this report. Thus, when 
the form of f, given by equation (8) is inserted into equation 
(5) for the local Mach number distribution and evaluated at 
the surface of the sinusoidal wall (y—0), the following result 
is obtained: 


M?—1 = H 3 
i1—AM 37 Tilt k 25 (n—2m)cos (n—2my » ; A* ^ 2 
© n=l mm - p=0 


(52) 
where H denotes the integral part of n/2. At the crest of 


the wavy wall (z--0), the point of maximum fluid velocity, 
with the numerical values for the coefficients inserted, equa- 
tion (52) becomes: 


w- 1 f ., 337 4043 
a) bg k +397 i+ oot ^ gz6x256 ^ V 
359381 ,,, 7326757 ,, 
270x250 © 6480x256 ^ + 
81688733 
86400x256 ^ T: (53) 


The critical value of & (that is, when M=1) calculated from 
this equation is 

k,,— 0.83244 l 
(Note that in ref. 1 the value k,,=0.83770 corresponds to the 
first six terms of eq. (53).) Consider now the infinite series 


= 1.60172 


4/3 
N * De3044 


LO 2.0 3.0 4.0 50. 6.0 ~ 7.0 
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Figure 2.—Numerical test of convergence. 
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i n=l 
where : 
a 3N— 9 
Aj n5) (n=2,3,... œ) 
ED 
The Cauchy ratio test A, — Au yields in the limit n— o 
i 0 0 


the result that the radius of convergence A, is equal to or 
less than 4/3. If the corresponding ratio Rn are formed for 


the right-hand side of equation (53) and the quotient N, =i 
2n 
is calculated, the resulting sequence of numbers is as follows: 


The noteworthy feature of this table is that, although 
Rin a presumably Kın) is converging quite slowly toward 


B (and R), the quotient N, exhibits a strong tendency 


to E an asymptotic value N for a relatively low value 
of n. Figure 2 shows this tendency in & graphic manner. 
The apparent asymptote represented by the straight line is 
the ratio of 4/3, the limit of R,, as n—>o, and of 0.83244, the 
critical value of b. Certainly, the rapid approach of the 
lower dotted curve toward the apparent asymptote and the 
decreasing oscillations represented by the upper dotted curve 
indicate that the critical value of k is the radius of conver- 
gence Rs of the power series on the right-hand side of equation 
(53). -Thus, the critical value of the stream Mach number 
marks the limit of convergence of the smooth symmetrical 
type of potential flow assumed. The ability to approximate 
closely the limiting value N is a matter of luck; namely, the 
choice of the known comparison series. Once, however, 
the proper comparison series has been selected and the 
approach to an asymptote indicated, there can be no question 
of the meaningfulness of the approximate value of N obtained. 
It may be that one would like to extend figure 2 to n=9. 
This extension would entail the forbidding calculation of an 
additional 185 coefficients A? ?. The result presumably 
would be to decrease slightly the approximate critical value 
of k and thereby raise slightly the straight-line asymptote of 
figure 2. This extension of figure 2 would show still more 
convincingly the approach to an asymptote and the dying-out 
of the oscillations. Perhaps more important still, figure 2 
definitely shows that conclusions based on less than six or 
eight terms are mere speculations in this field. 


CONCLUDING .REMARKS 


If the numerical test of convergence presented is accept- 
able, the conclusion to be drawn is that smooth symmetrical 
potential flow past the wavy wall exists only for the purely 
subsonic range. Moreover no such flow can possibly repre- 
sent the transonic or mixed type for which a local region of 
supersonic flow near the solid boundary is imbedded in the 
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otherwise subsonic stream. It follows as a corollary that 
the transonic or mixed type of flow past the wavy wall is 
necessarily an asymmetric one. This asymmetry in the 
flow pattern entails a resistance usually defined as wave 
drag. As shown by experimental observations, the shock 
wave associated with wave drag closes the downstream 
portion of the local supersonic zone. 

As a final remark—although the analysis and conclusions 


of the present work refer directly to the wavy wall, the’ 


suggested result that the critical stream Mach number 
marks the limit of smooth potential flow very likely applies 
to other boundaries. This conclusion is based on the idea 
that the gradual transition from a purely sinusoidal wall to 
& boundary composed of & single bump, say, introduces no 
essential changes in the analysis. 
LANGLEY ÁERONAUTICAL LABORATORY, 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS, 
LawcnEY Fignp, Va., May 6, 1952. 


APPENDIX 
`” NUMERICAL VALUES OF THE COEFFICIENTS 4*7 


Coefficients for Si: 
10... i13. 
Atom Aji N 256 
TEN d nı 2765 
451—556 27 9X 256? 
to. 1801 11 184345 
2 3X 256? 1728256? 
4896755 1 
10... IT. lll 
453—915 2563 dii 64 
33 
Il 
Ala= 64X256 
139 
11.— 
diz 4X256? 
Coefficients for fs: 
| la 2i — 125 
~ 10 As i= 36X256 
419 45895 
20— 2 1a a 
4i 72X256 As a 108X256? 
20234215 Aa 4210409 
277 432X256? $ 81X256? 
42791533 o 
20. alc 
As S= 648520568 dii 256 
20. l 21.815 - 
Ai =g di a= 576X256 
ao l 21 456185 
Ai V 1024 A S ^ 8450x250! 
4085 1 
20... ip —— 
Ai 2 7924x250? 4a i= 128 
21287 11 
20. Se 21... . 
Ay 3°" 324256? Aaa 32256 
Ario E _ 
ids 192x256? 


A= 65 13. 3386 
0 576X256 $ 1152x250? 
26435 , 3907 
Tg. o c ae tisdale 
Ao 52568 Ais 432X256? 
l 7 23 
Ti e ah ees 
Ai i= 72X256 Ais "12256? 
po 53995 Er | 583 — 
ps 1728x256? on 54x 256? 
9 119 
12.—.. 139—. ———É 
da= 32X256 E 12256? 
2765 ; 1 
12 . 13— m e—a 
Ai s= 96X256? As s: 96X256 
1 
12 — 
43 a= 8X 256 
55 
12. 
lm 8X256? 
a2. 12245 s3_____ 17999 
453—114 2561 = 90X256? 
A2 35023583 yIM 72577 
3 1620x256? 13 5760X256? 
Appa S28) 42a 45409 _ 
13 72X256? 33 2160x250? 
A22 104063 . 423 3515 
13 90X256? 33 864X256? 
79 | 10229 
29. 28. ——ÓÉÉÉ 
d3 3— 3X256? Ais 864X 2562 
15, 167 33. 79 
Az = 4256? As = 8256? 
47 25 
223 — 23. — 
dii 192X256 Ae § 12X256? 
VIE 24199 
337 864X256? 
A1i——- 
z 12X256 
11 
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Coefficients for f: 


1765 475823 
3 0— Sla ell 3 2__ 
455—753 Ai I= — 35556) | AS $— TES CO58S 
23603 1351393 74353 
A3 0— 31-—.. 323. 
01727 X 2562 xL 2880 X 256? Aii 320 256? | 
43 9... 9843883 ERE 15 | Asa... 3083 
93 91X 2563 p 32256 227 24X256? 
_ 1 30365 . Al 
Ai 9751 Ai a — 78x256 43i= 3x256 
259 117 1785 
30— $i... —_— AE 
E 72X256 4ii 32X256 A= 40X256? 
149317 29545 25 
30. $1. seit exe EN 
Als 432X 256? As; 96 X 256? Aii 2x 256? 
1 1 
30—. 7 Silcoci- 
xL 32 dzi 256 
33 55 
3 0— 31... 
4i 32256 Azi 256? 
139 
30. 
A11—3 58 
Coefficients for f,: 
7 35251 741115 
a EE 4d cc copo dg. Eo EE. 
990763256 Asi 90256? A53 5184 256? 
VIE 390547 . Asim 80872738 yIS 38381086 
720X256 675 X 256° 7405 2563. 
Aton 1704729613 pou 2465 ym 423833 
937 10200x 256? 11 " 576X256 22 1152X256? 
7 30596538 20845 
4$0—. ". p ag Oc erdndia Ded siad io CP MB ibo rsen iR 
Ai =g Aiz 135X256" © 433— 14 x 2568 
32947 455 173 
40 _ —  —— Ld ee E ORCSRN 42 . 
Aii 36 X 256? Aii 96X256 Aii 12256? 
20527213 127825 30 
40 L4 M 2 — d 
Aii 162 X 2568 Aii 192X256? Aii 2567 
3 l ^ 119 36 
t 0 41— .— CERES Ree 
Ais 128 Azi 48X256 Ags 5X256? 
4o, 017 At 1 ____ 91343 _ 
217 288 X 256 32 216x250? 
Akon 369587 Hite. 25 
23771728 K 256? = 48X256 
275 
T: At a= — 83359! 
1i 
40. 77 
Ait 24X256 
VT 4811 
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Coefficients for fs: ` Coefficients for fz: 
7 5470394 
Ag EET Y: Aj i= —2T5ER 70 13 
3840 81X256 Ao 0136X258 
90354659 88735 
Al SO Af i= ro d 
359200 X256 96x 256 191830286 
?9 144X256 
886243 25X997 l 
Ai t= Al = —5 359i jin oos 
1080x256 24x 256 Age 38852956 
25 105 
AS o= AS i= .. ~470_. 343 
1536 64x 256 44 — Tea O56 
24375 9 
2d. ; : Aj i= — e 70° 343 
32X256 32 X256 : Aj 0—320 x 286 
É - 
Åi o= agi E to. 2401 
di As 0435256! 
so... 8950 
31^ orny ore? 
27 X256* Coefficients for fj: 
AS Sm 25 : l 
to 24X256 482m 13X19 
° 18X256? 
VII 1375 " 
51 24X256? Aso 13X19 
, !9" 576X256 
Coefficients for fs: 
. A8 o. 89 — 
AS 0— 91 f PL 2483020 20 64X256 
° 432X256 t 54X256? l 
285522407 4371541 As I DUE 
60. 61... . A E 24X25 
E 48607256? Aii 5760 X 256? 
47 
91 265901 ETITI 
8 0— 6 1i ee 
Aio 12X256 E 192X256? xard 
i 19 
144862403 403535 A89-— 
Bo. eee eee. Eds eee EA 
A= 37 05¢2568 As 1— 388 x 256? .. 2880 
23 79841 j TRE 
60... S1. A8 0— —— 
39 83256 44 1—— 963x356: °° 2880 
2553091 2205 1 
60.— 61.. . S0 | 
dzi 2880 X 256? Asi 8X 256? Ar o= F599 
13 2401 
8 0— 61... — 
Aso 1024 E 60X256? 
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